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Introduction 


This report summarizes the work done in the period February 1993 through July 1993 on 
the “Prediction of Thermal Cycling Induced Cracking In Polymer Matrix Composites” 
program. An oral presentation of this work was given to Langley personnel in September 
of 1993. As per an oral agreement with Steve Tompkins, this presentation was accepted in 
lieu of the usual semi-annual report. This document was prepared for archival purposes. 

Progress 

Studies have been performed on the effects of spatial variations in material strength. 
Qualitative agreement was found with observed patterns of crack distribution. These 
results were presented to NASA Langley personnel in November 1992. 

The analytical methodology developed by Prof. McManus in the summer of 1992 (under an 
ASEE fellowship) has been generalized. A method for predicting matrix cracking due to 
decreasing temperatures and/or thermal cycling in all plies of an arbitrary laminate has been 
implemented as a computer code. The code also predicts changes in properties due to the 
cracking. 

Experimental progressive cracking studies on a variety of laminates were carried out by 
graduate student Cecelia Park at Langley Research Center. Results were correlated to 
predictions using the new methods. Results were initially mixed. This motivated an 
exploration of the configuration of cracks within laminates. 

A crack configuration study was carried out by cutting and/or sanding specimens in order 
to examine the distribution of cracks within the specimens. These investigations were 
supplemented by dye-penetrant enhanced X-ray photographs. The behavior of thin plies 
was found to be different from the behavior of thicker plies (or ply groups) on which 
existing theories are based. Significant edge effects were also noted, which caused the 
traditional metric of microcracking (count of cracks on a polished edge) to be very 
inaccurate in some cases. With edge and configuration taken into account, rough 
agreement with predictions was achieved. 

All results to date were reviewed with NASA Langley personnel in September 1993. View 
graphs from this presentation are attached. A paper on the previous work was also 
completed during this period and is attached. 1 




Current Status 


Cecelia Park is preparing her Master's thesis, which will fully document all of the work 
reported here. Work on a small thermal cycling and aging chamber is also in progress. 
The modified computer code is being documented and prepared for distribution. Work is 
commencing on the development of a modified theory to account for crack configuration, 
strength variation, and edge effects. The paper 1 , recendy presented at the American 
Society for Composites 8th Technical Conference on Composite Materials, is being 
prepared for submission to an archival journal. 
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ABSTRACT 

Thermal fatigue has been observed to cause matrix cracking in laminated composite 
materials. A method is presented to predict transverse matrix cracks in a composite laminate 
subjected to cyclic thermal load. Shear lag stress approximations and a simple energy-based 
fracture criteria are used to predict crack density as a function of temperature. Prediction of 
crack density as a function of thermal cycling is accomplished by assumin g that fatigue 
degrades the material’s inherent resistance to cracking. The method is implemented as a 
computer program. Simple experiments provide data on progressive cracking of a laminate 
with decreasing temperature, and on cracking induced by thermal cycling. Correlations of the 
analytical predictions to the data is very good. A parametric study using the analytical method 
is presented which provides insight into material behavior under cyclical thermal loads. 

BACKGROUND 

The orbital environment is characterized by wide swings in temperature as the vehicle 
moves in and out of the shadow of the earth. Composites subjected to this thermal cycling 
suffer from transverse matrix cracks [1]. These cracks cause degradation of material 
properties, dimensional changes, and changes in the coefficients of thermal expansion (CTE) 
of the composites. In extreme cases, this can cause the material to simply come apart More 
relevant is the changes in properties that can be caused by even modest levels of cracking. 
Composite laminates are often used in space structures due to their extreme dimensional 
stability. Such laminates are, however, very sensitive to imperfections, and their CTE's can 
be changed significantly by matrix cracking levels that are not otherwise a threat to the 
structure. This problem is under continuing experimental study. 

The problem of transverse matrix cracking in laminated composites due to mechanical 
loading has been extensively studied. It has been observed that Classical Laminated Plate 
Theory (CLPT) cannot predict transverse cracking; the strain at which transverse cracks 
appear is dependent on the thickness of the cracking ply-group [2], a dependency which is not 
predicted by CLPT. In practice, the transverse strength of a group of plies with the same 
orientation is often treated as a function of the ply group thickness. To be used in design, 
however, this approach requires data from many different layups for each material under 
consideration. A predictive methodology based on a shear lag model of the stress distribution 
around a transverse crack, combined with an energy-based criteria for crack growth, has been 
studied by many authors. A good review, which will not be repeated here, is given by Naim 
[3], and extensive work continues to appear in the literature [4,5]. The method of Laws and 
Dvorak [6] is particularly well suited to investigations of progressive damage due to 
monotonically increasing loads. 
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Transverse cracking has been recognized as the first stage of damage due to cyclic loading 
in many laminates [7-9]. The literature concerned with the characterizadon and prediction of 
fatigue damage has extensive descriptions and measurements of this phenomena, but is less 
concerned with predicting it than using it as a metric of damage. In particular, reference is 
made to a characteristic state" of damage at which the laminate is allegedly saturated with 
cracks. The crack distribution at this state is taken to be a laminate property, independent of 
loading type or history. 

A limited literature exists on attempts to use shear lag and fracture mechanics concepts to 
predict transverse cracking as a function of fatigue loading. Petitpas et al. [10] use a shear lag 
stress solution and a maximum stress criteria for crack formation, coupled with measurements 
showing that the effective transverse failure stress is a function of number of cycles, to predict 
the crack density as a function of cyclical loading. Lafarie-Frenot and Henaff-Gardin [11] 
propose an empirical relation between crack density and number of cycles in the form of a 
Paris-type law. 

Most of these studies include the important effects of residual curing stresses, which are 
thermal stresses. None, however, explicitly address the issue of transverse crack formation 
under thermal loading. 

STATEMENT OF PROBLEM 

Given appropriate material properties, the layup and geometry of a composite laminate, 
and a thermal loading history, we predict the resulting distribution of matrix cracks, and the 
resulting degraded laminate properties. An analytical model is implemented as a computer 
program to provide the predictions. It is verified by correlation with experiments, some 
performed as part of this effort and some garnered from the literature. 

ANALYTICAL MODEL 

Thermal or mechanical loading stores strain energy in the laminate. Some of this energy is 
released by the formation of a transverse crack. Typically, these cracks appear suddenly 
across the full height and width of a ply or ply-group (group of adjacent plies all with the same 
alignment angle) parallel to the fibers. Figure 1. The amount of energy released can be 
approximated by a shear lag solution of the stresses in the vicinity of the crack. The 
appearance of a crack can then be predicted by a simple energy method- if the energy released 
AG is greater than the energy necessary to form the crack G fc , the crack will form. 

A simple shear lag model was derived. It was found to be identical in form to that of 
Laws and Dvorak. Readers are therefore referred to Reference 6 for details. We will use the 
notation of Laws and Dvorak with minor variations. We consider a laminate consisting of a 
ply group with fibers aligned 90° to the x axis embedded in an arbitrary laminate. Figure 2 
shows an edge-on view of such a laminate. A unit depth in the y direction is assumed. We 
assume the cracks labeled A and B exist, and wish to find the strain energy released by the 
appearance of the additional crack C. The thickness of the cracking ply group is a 2 \ a, is the 
combined thickness of the remainder of the laminate, and a 0 is the laminate thickness. 
Similarly, E 2 is the uncracked axial stiffness of the cracking ply group, E, is the axial stiffness 
of the remainder of the laminate, and E„ is the stiffness of the entire laminate. These effective 
stiffnesses can be calculated from the ply properties and laminate geometries using the 
relations contained in the Appendix. This is, with minor differences in notation, the problem 
studied by Laws and Dvorak. They derived 

a,{oy+cu)E g L £/t, L ^2^] ,,, 

2£a,£,£, V <h <h <h y 

where AG is the change in strain energy due to the appearance of crack C. | is a geometric 
shear lag parameter, and a is the axial stress in the cracking ply group prior to cracking. 

Laws and Dvorak expressed this stress as a function of applied load and residual thermal 
stresses. We wish to express it in terms of applied thermal load AT. In the absence of applied 




Figure 1 . Cracked cross-ply laminate 


Figure 2. Cracking geometry 


mechanical load, the stress due to a change in temperature can be shown to be 

r _ a l 5&te-oQA7’ 


( 2 ) 


where AT - T - T !f and 7> is the stress-free temperature of the laminate, often assumed to be 
the cure temperature. In materials of interest, a 2 > a h so a negative AT is required for a 
positive stress. Combming Eqs. 1 and 2, and setting AG to G lc , the temperature change 
necessary to form a crack at location C is found to be 


Ar = _ 

a i a 2 E l E l(**7-<Xl) 2 t 


tanh^ + tanh-^-tanh^^- 
V <h <h 


(3) 


The negative value of the square root is selected on the physical grounds that this value of AT 
produces positive stress in the cracking ply group. 

Multiple cracks are handled by predicting a crack density p, which is the inverse of the 
average crack spacing. Predictions based on a single parameter p assume an unrealistically 
even crack spacing, but can be used to calculate either the maximum possible crack density, or 
e minimum density. For example, if we assume a uniform crack spacing, and assume new 
cracks will form half way between the existing cracks (h,=h 2 =h in Figure 2), Eq. 3 becomes 


AT* = _ + ) 

i ( h<L l E x E 2 {a 2 -a x f t 


2tanh— - tanh-^^- 


*2 


Oy 


(4) 


An existing crack spacing of exactly 2h will allow many more cracks to form at AT*, 
resulting in crack spacing h; while an existing crack spacing of just under 2h will not allow 
an>r m°re cracks. Therefore, the minimum possible crack spacing due to temperature change 
ni is h, while the maximum possible spacing is just under 2h. 

The above formulation does not account for the observed effects of thermal cycling To 
account for these effects, we modify the approach of Petitpas et al. [10]. They report a 
decrease in the first ply failure stress of cross-ply laminates due to mechanical loading cycles. 
By modifying Eq. 1, with h - h t - h 7 ~ <», and setting AG to G, e , we can relate the stress in the 
cracking ply group at the time of the first crack <3r c to the cridcal strain energy release rate. 


2ME, 




( 5 ) 


Given effective failure stress as a function of loading cycles N. and the original critical 
stain energy release rate, we can reasonably estimate the effect of cyclical loading on the 




critical strain energy release rate by dividing Eq. 5 by the zero-cycle values and rearranging. 
G, c (W) = C„(0)[<T Jt (W)/<T /c (0)] ! ( 


Given G, c as a function of number of loading cycles, we can predict the limits of crack 
density as a function of constant thermal cycles by solving Eq. 4 for A, with AT* the absolute 
maximum AT during the cycle, and G lc being G k (N). This procedure ignores stress ratio 
effects. They could be included without change in the procedure by using aJN) data 
collected using the correct stress ratio. The procedure also presumes that all material 
properties are temperature independent 

The degraded properties of the cracked laminates are now calculated as functions of crack 
density. Laws and Dvorak [6] give a general derivation of this problem. We will start with 
their formula for reduced laminate axial stiffness, £* . Using our notation, and noting p = 1/h 


£r=£„/ll + ^^-tanh^- 

2^ pa, J 


(7) 


Clearly, the reduction in laminate stiffness is due to a reduction in the stiffness of the 
cracked ply group. If E'_ is this reduced stiffness 

= kE, (g) 

where vis a knock-down factor due to the cracks. It can be calculated by combining Eqs. 7, 

8 and A4 » • 


v=£;a, 


2£ 


P<hJ/ 1 


pc u 

24 



(9) 


All of the properties of the cracked laminate can now be calculated by CLPT [12]. The 
properties of the cracked plies are modified by multiplying the non-axial reduced stiffnesses 
by the factor jc The details of these calculations are found in the Appendix. 

COMPUTER PROGRAM 


The models of laminate behavior outlined in the previous section were implemented as a 
computer code, CRACK-O-MATIC. The code takes material properties, laminate geometries, 
and thermal loading histories, and predicts minimum or maximum crack densities, and the 
corresponding degraded laminate properties. The code currently performs three types of 
calculations- 

1) prediction of the crack density due to monotonic application of a given AT 

2) tabulation of crack densities and degraded material properties as a function of 
monotonically decreasing temperature 

3) tabulation of crack densities and degraded material properties as a function of 
number of constant thermal cycles 

The code is available in user-friendly form, including instructions, from the authors. 
EXPERIMENTS 


Two types of experiments were earned out as part of this effort. A group of specimens 
were exposed to steadily decreasing temperatures, and several groups of specimens were 
subjected to constant thermal cycles. The decreasing temperature specimens were made from 
P75/934 material. The thermal cycling laminates were made from prepreg containing 
continuous graphite fibers of different types and the ERL 1962 epoxy resin from Amoco 
Performance Products, Inc. All of the composite laminates were fabricated at the NASA 
Langley Research Center using vendor recommended procedures. Specimens were cut from 
12 inch square laminated panels of each materiaL One side of each specimen was polished so 
damage could be observed. All specimens were dried to a constant weight in a warm vacuum 



oven and stored in a vacuum at room temperature to eliminate moisture absorption effects. 

The layups and materials used are summarized in Table 1. All specimens measured 3.0 inches 
by 0.5 inches. Three replicants were used in each test 

The first test measured crack density as a function of decreasing temperature. The 
specimens were cooled to a target temperature at no more than 10°F/min, held at that 
temperature for 15 minutes, and returned to room temperature at no more that 10°FAmn. A 
small chamber equipped with liquid nitrogen cooling and electrical resistance heating was 
used. The specimens were shielded from contact with the heating and cooling elements, and 
were heated and cooled by circulating air only. After each excursion, crack densities were 
measured at room temperature by counting cracks under optical magnifications from lOOx to 
400x. The number of cracks in the middle one inch of each specimens was counted, and the 
results averaged. Only cracks that were observed to go more that half way through the central 
90 ply group were counted. Measurements were taken of the specimens as-received, and 
after excursions to 0°, -50°, - 100°, - 150°, -200°, and -250°F. The results are shown in the next 
section. 

Thermal cycling test were then performed. The damage state of each specimen was 
determined prior to thermal cycling by examining the polished edges with an optical 
microscope at a magnification of 400x. Specimens were thermally cycled between +/-250°F in 
a chamber similar to the environmental chamber described in Reference 13. For these tests 
the thermal cycle had a period of about 25 minutes. After specimens had been exposed to 
predetermined number of cycles, the specimen edges were examined for induced damage and 
then subjected to additional thermal cycles. Cumulative damage was tracked for each 
specimen. Again, the results are in the next section. 

RESULTS AND CORRELATIONS 

The results of the monotonic cooling experiments are shown in Figure 3. The measured 
crack density is shown as a function of temperature. The triangles represent the means of 
three specimens, while the bars represent the limits of the individual measurements. It can be 
seen that the data is very consistent. The experiment was not a true monotonic cooling test, as 
the specimens were returned to room temperature for crack counting, but the effect of these 
few thermal cycles were assumed to be negligible compared to the damage done by exposure 
to successively lower temperatures. The fact that cracks were observed in as-received 
specimens indicates that the cooling from the cure temperature of 350°F (which was assumed 
to be the stress free temperature) to room temperature was sufficient to initiate cracking. 

The code was used to generate predictions of laminate cracking. The ply property data 
used is in Table 2. The material stiffness properties are NASA Langley data. The strain 
energy release rate G/ c and the shear lag parameter £ for the P75/934 material was not known, 
so the values were adjusted to fit the data. The correlation is shown in Figure 3. The match 
between the analysis and data is excellent This is, of course, not surprising given the 
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Figure 3. Monotonic cooling data and analysis 



Table 1. Tests Performed 


Test Type 

Material 

Layup 

Thermal 
Cycling 
(+/- 250°F) 

T50/ERL1962 

T50/ERL1962 

P5 5/ERL 1962 

P5 5/ERL 1962 

P75/ERL1962 

P75/ERL1962 

P100/ERL1962 

PI 00/ERL 1962 

P120/ERL1962 

P120/ERL1962 

[0/+-45/90)s 
[0/90/0/901 s 
[0/+-45/90]$ 
[0/90/0/90] s 
[0/+-45/90] s 
[0/90/0/90] s 
[0/+-45/90] s 
[0/90/0/901 s 
[0/+-45/90]s 
[0/90/0/901 <; 

Monotonic 

cooling 

P75/934 

(02/902ls 



Table 2. Ply Properties 


Material 

E, 

Ml 

HI 

mm 

Y, 

MSM 

WEm 

KM 



(Msi) 




fksi) 

mssm 

Wsumam 

■filial 


P75/934 

34.3 

0.9 

0.29 

0.7 


40 

-0.68 

16.0 

0.65 

T30Q/934 

- 

- 

- 

• 


250 

. 


0.90 

T5Q/ERL1962 

28.8 

1.04 

0.27 



132 

-0.30 

18.0 

0.90 

P55/ERL1962 


1.00 

0.34 


3.25 

73 

-0.38 

16.0 

0.65 

P75/ERL1962 

34.3 

0.96 

0.29 

0.70 

3.88 

104 

-0.53 

■ 

0.65 

P10Q/ERL1962 

51.7 

0.87 

0.30 

0.70 

2.67 

50 

-0.53 

KISH 

0.65 

P120/ERL1962 

58.5 

0.87 

0.28 

0.70 

2.89 

58 

-0.68 

16.0 

0.65 


adjustments of the constants, but the constants cannot change the overall trend predicted by the 
analysis, so the fact a fit was achievable is at least encouraging. 

More definite verification of the utility of the analysis can be seen in the fatigue study 
correlation. In this study, two different layups and five different materials were studied. The 
analysis used material properties measured at NASA Langley and shown in Table 2. G, c 
values were not available, so they were scaled from available transverse ply strength dam, and 
the available G, e of T300/934 material. Assuming transverse strength Y, is a fracture 
dominated property, it is reasonable to assume that for mate rials r and s, 

( 10 ) 

It was found that this method produced values within the range of measured G, e values 
available for some materials. The value of the shear lag parameter was held constant for all 
fibers of the same type, although it was found by comparison to available daxa on T300 
materials that a value of 0.9 was appropriate for PAN based fibers, while the above correlation 
suggests 0.65 for pitch based fibers. Given the geometric differences between the fiber types, 
and the fact that the shear lag parameter is a geometric parameter, this is reasonable. 

The dependencies of G, c on the number of loading cycles for all materials considered were 
scaled using Eq. 6 and the data of Petitpas [ 10]. The result for all materials considered was 


G te {N) = G, c (0) (1 -0.0336] 2Io *io n 


(ID 


The results are shown in Figure 4. Crack density is plotted vs. number of cycles, with the 
mean of the data shown as symbols and the predictions as lines. The correlations are 




Cycles (+/-250°F) 


* 


50 


Layup 

Data Prediction 

[0/+-45/90 ]s 

A 

[0/90/0/90] s 

o 


Figure 4. Thermal fatigue data and analysis for several materials and layups 





















startlingly good. No modifications of the any of the analytical inputs were made to achieve 
these correlations. It can be seen that the effects of variations in material types and layups are 
successfully predicted, both qualitatively and, with the sole exception of the P55 material, 
quantitatively as well. Most notable is the difference between the two different layups of T50 
material- this subtle difference causes one laminate to crack after a few cycles, while the other 
remains uncracked, and the effect is successfully predicted by the code. 

PARAMETRIC STUDIES 

Once a reasonable level of confidence in the analysis was achieved, an attempt was made 
to obtain greater physical understanding of the cracking problem through parametric exercises 
with the code. Simulations of monotonic cooling and thermal cycling were performed, 
varying material critical energy release rate, cracking layer thickness, and surrounding 
laminate layup. The effect of crack densities on the laminate properties were also examined as 
a function of cracking layer thickness. 

Figure 5 shows predicted cracking due to cooling for a family of layups with different 
assumed G lir A [0V90 : ], layup and P7 5/934 elastic properties used for the calculations. 
Increasing G /c decreases both the temperature at which cracking initiates and the level of 
cracking at any given temperature. Note that in no case does the level of cracking reach a 
steady state- at least in theory there is not "characteristic damage state" for this laminate 
suffering damage due to monotonic cooling. 

Figure 6 shows the effect of decreasing cracking layer thickness. The baseline [0 2 /9O0 s 
laminate, with four 90° plies in the cracking group, is compared to a [02/90], (with 2) and a 
[02/90],, where the bar indicates that the laminate is symmetric about a single 90° ply. Thinner 
ply groups delay the onset of cracking, but once cracking starts it continues to higher densities 
in the thinner ply groups. Figure 7 shows the same study applied to a thermal cycling case. 
Due to the severe +/-250°F thermal cycles, the thicker ply groups are severely cracked on the 
first cycle. The laminate with a single (5 mil) center ply shows somewhat different behavior; 
if the single ply is reduced only slightly in thickness, to 4 mils, the laminate does not initially 
crack, but does crack after thermal cycling. This is in qualitative agreement with observed 
behavior [14]. 

Figure 8 shows the effect of different surrounding laminates on the cracking of a two-ply 
central 90° ply group in thermal cycling. The effect is noticeable, and the trends non-intuitive. 
The number of zero degree plies in the [0^/90], family of laminates (where n from 1 to 3 was 
checked) makes little difference. The more compliant [0/+-45/90], laminate has worse 
cracking, but the presence ot other 90° ply groups in the [0/90/0/90], laminate seems to help. 
Again, these results are partially confirmed by experiment- note the T50 results in Figure 4. 

The effects of cracking on laminate stiffness and CTE are seen in Figures 9 and 10. 

Figure 9 shows normalized stiffness loss. The stiffness loss is always fairly small, and is 
negligible in the case of thin cracking ply groups. Figure 10 shows the laminate CTE; 
normalization would be confusing in this case as the uncracked CTE's of each laminate are 
noticeably different Cracking causes the CTE to change more than 100% of the nominal 
value, and even change signs, in some cases. Obviously this is an important effect especially 
in structures requiring extreme dimensional stability. The change in CTE due to cracking is 
larger than the expected variation due to manufacturing and material variations [15], and hence 
may cause the structure to drift outside of design limits for thermal deformations. This effect 
is much less severe for thin cracking layers. 

CONCLUSIONS AND RECOMMENDATIONS 

The method of shear lag analysis combined with an energy method for crack prediction can be 
used to predict thermal cycling induced matrix cracking in composite laminates. The method 
has been verified by comparisons with several different kinds of experimental results. A 
parametric study using the analytical method gives good insight into the physical problem. 

One important result is that laminates with thinner ply groups crack at lower temperatures 
and/or larger number of thermal cycles, but when they do crack they will eventually reach 
higher crack densities. On the other hand, the effect on laminate properties when thin ply 
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Figure 5. Sensitivity of crack density to 
G k in monotonic cooling 
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Figure 7. Sensitivity of crack density to 
cracking layer thickness in thermal cycling 



Crack Density (per in) 

Figure 9. Stiffness loss due to cracking 
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Figure 6. Sensitivity of crack density to 
cracking layer thickness in cooling 
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Figure 8. Sensitivity of crack density to 
surrounding layup in thermal cycling 
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Figure 10. CTE change due to cracking 






groups crack is much smaller that the effect caused by cracking of thick ply groups. It was 
also found that the material's transverse strain energy release rate is important to the onset and 
eventual density of cracking, and that higher values are desirable. Less intuitive was the 
finding that the layup has an important effect on the level of cracking, but that no one factor 
dominates this effect 

These results are significant in light of recent efforts to reduce matrix cracking in 
composite space structures through the use of thinner prepreg plies [14]. The hope has been 
to eliminated cracking entirely. Our results indicate that even if no cracks are evident after a 
few thermal cycles, cracks may form after further thermal cycling, and ultimately may reach a 
higher density that would be reached with thick plies (see Figures 6 and 7). However, these 
cracks have little effect on the laminate properties, as seen in Figures 9 and 10 Thin plies 
may be useful for enhancing the stability of space structures, but for different reasons than 
those originally proposed. In is also evident that, when comparing laminates with different 
ply group thicknesses, crack density is alone is not a useful damage metric. 

Further work is needed to make this method useful for quantitative, predictive analyses 
Accurate critical strain energy release rate data, preferably over a wide range of temperatures 
and collected using a variety of specimen types and materials, is required. The method also 
needs to be generalized so that cracking in more than one layer, and in layers with ply angles 
of other than 90°, can be predicted. Effects such as temperature dependence of material 
properties, stress ratio effects, and random load cycles also need to be included. 
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APPENDIX 

We assume that the cracked laminate is much wider than it is thick. Therefore, the 2- 
dimensional problem presented in Figure 2 is a plane strain problem. The appropriate 

m^^Hni UffneSS !f- E °'r h ,^u h ** calculated hy CLPT [12]. For each ply t, we have 
measured ply properties £**, (fiber direction stiffness), (transverse stiffness) v (maior 

Poisson s ratio), Cj (shear stiffness), (fiber direction CTE), and a* (transverse CTE). The 

ply has thickness r„ and the fibers are aligned at an angle ft to the x axis. The cracking dIv 

gr ? U ^ ^ a sm J e P‘y k ' wherc ft = 90°. The equivalent stiffnesses are then 

calculated as follows. For each ply i 
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For the laminate 
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and we recalculate Eqs. A1-A3 using these values in Eq. Al when assembling Q k . 
effective stiffness E« and effective CTE cr^ofa laminate are defined as 

E tff = 1 

A' = A" 1 
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• treater understanding or this problem is necessary. Want 
to be able to predict matrix cracking and its affect on 
laminate properties 



BACKGROUND 


Crack Formation In 
Thermally-Loaded Composites 




Constrained Laminate 


Cracks develop due to 
thermal (CTE) mismatch 





CLOSEUP VIEW OF MICROCRACKING BEHAVIOR 




PREVIOUS ANALYTICAL WORK 
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Our problem requires thermal, progressive and fatigue 
damage study for all types of laminates 



SUMMARY OF WORK BY H.M. 
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• Overall, good results for analysis of crossplies 



RESEARCH OBJECTIVE 
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For all plies and arbitrary layup 



GENERAL TASKS: 
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thermal progressive and cycling experiments 



EXPERIMENTAL PROCEDURE 

Progressive Thermal Loading 

Material: P75/934 Graphite Epoxy 
Size: 0.5 X 2.5 inches P y 

Layups: [02/902]s Quantity: 3 each 

[0/+45/90/-45] s 9 total 

[0 2 /±30] s 

t0 P^^ssively lower temps. 
( to 200 F) in thermal environmental chamber 

£^ ac ^ d u en i it ^ (oracks/inch) and distribution taken 
on polished edge (SideA) using optical microscop! 



Time (not to scale) 




EXPERIMENTAL RESULTS 
PROGRESSIVE LOADING 
[02/±30]s 



Temperature (°F) 


• Specimens hardly cracked, some cracking 
in +30 plies 

• Preliminary analysis predicts zero cracks 




EXPERIMENTAL RESULTS 
PROGRESSIVE LOADING 
[0/+45/90/45] s 



Temperature (°F) 


• Trends for the middle -45 plies match the 
preliminary analysis 

• Preliminary predictions for 90° plies start to 
crack too late but end up with same density 

• +45° plies start at correct temperature but 
crack density is low 

• 90 and +45°data does not seem to agree 



CRACK CONFIGURATION STUDY 
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2) Inspection of other edge (Side B). 

3) Sanding edges 



XRAY PHOTOGRAPHY 


RESULTS: 


[02/902ls : Cracks propagate through entire width of 
specimen. About one crack/inch does not go through. 



[0/+45/90/*45]s: Thick middle -45 layer showed good 

behavior like the crossplies but pictures not as clear. 
Evidence of some cracks in other plies after cooling 

[0 2 /±30] s : Inconclusive 


CONCLUSIONS: 

• Cracks in [0 2 /90 2 ] s specimens behave 'ideally' 

• Except for middle -45 layer, cracks either too thin/small 
to see on photographs or do not propagate through the 
width 

• [0/+45/90/-45] s and [0 2 /±30] s specimens need further 
examination 


EDGE COMPARISON OF SIDE A AND B 



• Some cracks in thin plies do not propagate through. 

• Cracking behavior is dependent ply thickness 




SANDING EXPERIMENT 
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• Crack density and distribution taken for [0 2 /±30] s and 
[0/+45/90/-45] s specimens 



SANDING EXPERIMENT RESULTS 
[0 /±30]s SPECIMENS 
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CRACK DENSITY (PER INCH) CRACK DENSITY (PER INCH) 


SANDING EXPERIMENT RESULTS 
[0/+45/90/-45] s 
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+ 45° AND -45° PLY DATA 



0 2 4 6 8 10 12 


90° PLY DATA 



0 2 4 6 8 10 12 




procedures used so far (fine, rough and fine/rough combo) 

Cracks will only appear in new 'undamaged' specimens if 
sanding causes the damage. 



CONTROL TESTS 
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SANDING RESULTS USING FINE SANDPAPER 

[0/+45/90/-45L 
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DISTANCE ALONG WIDTH OF SPECIMEN (mm) 




SANDING EXPERIMENT 
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OVERALL CONCLUSIONS FROM 
CRACK CONFIGURATION STUDY 



• Must be cautious in using edge crack density count as 
measure of damage 
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